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Abstract 

We prove a recently conjectured star-star relation, which plays the role of an integrability 
condition for a class of 2D Ising-type models with multicomponent continuous spin variables. 
Namely, we reduce this relation to an identity for elliptic gamma functions, previously obtained 
by Rains. 
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1 Introduction 



Recently two of us [Tl[2] introduced a new class of exactly solvable 2D lattice models of statistical 
mechanics, which involve continuous spin variables taking values on a circle. The interest in these 
models is motivated by various applications. In statistical mechanics they serve as rather general 
"master models", containing many important particular limits, such as the Ising, chiral Potts Odl, 
Kashiwara-Miwa [5], Faddeev-Volkov [6l|7] and other models. Mathematically, the new models are 
deeply related to the theory of elliptic hypergeometric functions [8]. For instance, the celebrated 
elliptic beta integral [9], which lies at the basis of this theory, is shown [1| to be a Yang-Baxter (star- 
triangle) relation, defining perfectly physical integrable lattice models of statistical mechanics. Other 
interesting connections are discussed in [10^13]. We mention, in particular, that the elliptic gamma- 
functions arise in calculations of superconformal indices connected with electric-magnetic dualities 
in 4D J\f = 1 superconformal Yang-Mills theories [14] . Most remarkably, as recently discovered 
in |15H17j , the superconformal indices in 4D superconformal quiver gauge theories precisely coincide 
with partition functions of the 2D lattice "master models" [HE] discussed here. Interestingly, 
in this correspondence the Seiberg duality for the superconformal indices reduces to Baxter's Z- 
invariance |18j for the partition function under (generalized) "star-triangular moves" of the 2D 
lattice. 

Here we resolve an outstanding question for these models concerning the so-called star-star 
relation, conjectured in [2]. This relation serves as an integrability condition, as it implies the 
Yang-Baxter equation, the Z-invariance of the partition function and the commutativity of row-to- 
row transfer matrices. In this letter we completely prove this star-star relation by reducing it to a 
transformation formula for elliptic hypergeometric integrals, previously obtained by Rains |19j. 



2 Edge-interaction model with continuous spins 

In this Section we formulate the star-star relation conjectured in [2] (see (I14p below). First, we need 
to briefly describe the associated two-dimensional edge-interaction model; full details can be found 
in [2]. Consider the regular square lattice, drawn diagonally as in Fig. [TJ The edges of the lattice 
are shown with bold lines and the sites are shown with either open or filled circles in a checkerboard 
order. At the moment we will not distinguish these two type of sites. At each lattice site place a 
n-component continuous spin variable 

n 

a; = {xi, . . . ,Xn} G M", < Xj < vr, '^Xj=^ (mod vr) . (1) 

i=i 

Note that due to the restriction on the total sum, there are only (n — 1) independent variables Xj. 
For further reference define the integration measure 

dx = ■ ■ ■ dxi ■ ■ ■ dxn-1 , (2) 
Jo Jo 

over the space of states of a single spin. Fig. [T] also shows an auxiliary medial lattice whose sites 
lie on the edges of the original square lattice. The medial lattice is drawn with alternating thin 
and dotted lines. The lines are directed as indicated by arrows. To each horizontal (vertical) line 
on the medial lattice assign a rapidity variable u (v). In general these variables may be different 
for different lines. However, a convenient level of generality that we shall use here is to assign the 
same rapidity u to all thin horizontal lines and the same variable u' to all dotted horizontal lines. 
Similarly, assign the variables v and v' to thin and dotted vertical lines as indicated in Fig. [TJ 
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Figure 1: The square lattice shown with bold sites and bold edges drawn diagonally. The 
associated medial lattice is drawn with thin and dotted horizontal and vertical lines. The lines are 
oriented and carry rapidity variables u, u' , v and v' . 



Two spins interact only if they are connected by an edge. To define the Boltzmann weights we 
need to introduce the elhptic gamma- function [8]. Let q, p be two elliptic nomes (they play the role 
of the temperature- like parameters), 

p = e''^'^ , q = e'''^, Imcr > 0, ImT>0, (3) 

and 

r/ = -i7r(cj + r)/2 , (4) 
denote the "crossing parameter" . Define the elliptic gamma-functior0 



- I_g2i.q2,+lp2fc+l 

j,k=o' ^ "I P 



^ Q—2\zk 



where the product formula is valid for all z, while the exponential formula is only valid in the strip 

— Re r/ < Im z < Re 7/ . (6) 

The function ([5]) possesses simple periodicity and "reflection" properties 

$(z + 7r) =$(z), <S>{z)^{-z) = I. (7) 

Each edge is assigned a Boltzmann weight which depends on spins at the ends of the edge and on 
two rapidities passing through the edge. There are four types of edges differing by orientations and 
types of the directed rapidity lines passing through the edge. They are assigned different Boltzmann 
weights as shown in Figl2j These weights are defined as 

n 

Wa{x,y) = Kn{a)-^ Yl ^{xj-yk + \a), Wa{x,y) = y^S{x)E{y)Wr,-a{x,y) , (8) 
j,k=i 



^Our function ${z) coincides with r(e ^''"^ p^, q^) in the notation of ref. [8]. 
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Figure 2: Four different types of edges and their Boltzmann weights. 



where Kn{a) is a normalization factor. The single-spin function § is given by 



§{x) = Jl |$(xj -Xfc + iry)} 



vr 



n-l 



G(q)G(p), 

where the indices j, k run over the values 1,2, ... ,n and 

oo 

G(z)=n(i-^'')- 

A;=l 

The partition function is defined as 

Z = j \{mu-v'{x,y) []w„_„(y,x) \{Mu'-v{xy) J] W,._,.(y, a;) \{dx. 



{xy) 



{y,x) 



sites 



(9) 



(10) 



where the four products are taken, respectively, over the four types of edges shown in Fig. [2j The 
integral is taken over all configurations of the spin variables on the internal lattice sites. The 
boundary spins are kept fixed. 

Note that the lattice in Fig. [T]can be formed by periodic translations of a four-edge star, con- 
sisting of four edges meeting at the same site. A little inspection shows that there are only two 
different types of such stars shown in Fig. [3l They are either centred around "white" sites, shown 
with open circles or around "black" sites, shown with filled circles. Applying the rules of Fig. [2] one 
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Figure 3: Two types of four-edge stars: a white-centred star V*-^' (left) and a black-centred star 

V(2) (right). 
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can write Boltzmann weights corresponding to these stars 



VL^2 " d ) = / dxWu-v{c,x)Wu'-v'{b,x)Wu'^v{x,a)Wu-v'{x,d) , (12) 

and 

^^'^ ( c d ) " / dy^u-v{y,b)Wu'-v'{y,c)W^,^,{d,y)Wu-v'{a,y) , (13) 

where the bold symbols u = [u, u'] and v = [v, v'] stand for the rapidity pairs. It turns out that the 
above two expressions are simply connected to each other 

W,,_,(d,c)W„,_Jd,6) vL'i " ^ ^ = W,._,(b,a)W,._Jc,a) " (14) 

This is precisely the star-star relation conjectured in [2] (in the same paper it was also verified in a 
few orders of perturbation theory in the parameters p and q). Below we will give a complete proof 
of (I14p by reducing it to a mathematical identity, previously obtained by Rains [l9j. Apparently, 
the above star-star relation is the simplest condition for the Boltzmann weights which ensures the 
integrability of the considered modell. In particular, it implies the commutativity of the row-to-row 
transfer matrices. 



3 Proof of the star-star relation 

Here we will use the standard notation for the elliptic gamma-function, which is simply related to 
our definition 

^(z;p^q2)=^ ^ ^.(x) = ^(pqe-2^p^q^), pq = e'^^ (15) 

j,k=o ^ p q 

In the following we will omit the nome arguments p^ and q^, assuming that r(2;) = r(2;;p^,q^). 
Following Rains [19] introduce the following elliptic hypergeometric integral! 

= ;^(^;^j — n,^,r(../.,) — Hir- 

which involve 4n independent parameters 

{ti} = {tl,t2, ■ ■ ■ ,t2n}, {Si} = {si,S2, . . . ,S2n}, \ti\,\Si\<l, i = I, . . . ,2n , (17) 

where n > 2 and the function G{z) is defined in (jlOp . The indices k and i in the denominator of 
(I16p run over the values 1, . . . ,n. All integrations are taken over the unit circles |zfc| = 1, and the 
variable z„ is determined by the constraint 

ZlZ2 ■■■Zn = l. (18) 

^For n = 2 the star-star relation (I14|l is just a consequence of the star-triangle relation, Eq.(1.5) of [T], which is 
equivalent to the elliptic beta integral 9, 10 . However, for n > 3 the corresponding star-triangle relation apparently 
does not exist (at least it is not known to the authors) and the star-star relation H14p seems to be the simplest relation 
of this type. 

■^We follow Sect. 4 of [19], where we set Z = 1. Moreover, our indices numerating parameters in ()17p start from 1 
instead of in [19j . 
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The theorem 4.1 of [T2] (where we set m = n and Z = 1) states the fohowing transformation 
formula, 

2n 

{s.}) = ( n ntj^k)) {s^}) (19) 

j,k=l 

where the new parameters {ti} and {sj} in the RHS are given by 

2n 2n 

ti = TT^ti-\ Ui = U^sr\ T = Yltj, U = Y[sj, i = 1,2, . . . ,27i . (20) 

i=i j=i 

Consider now the star weight (|12p and make a change of variables Zj = e+^'^J and 

f , _ p-2{u-v)-2\Cj . -2{u'-v')-2\bj _ 2{u'-v-T])+2\aj , _ 2{u-v' -ri)+2\dj 

where Oj, 6j, Cj,dj are the components of the spin variables and j = 1, . . . ,n. Note that due to ([T]) 
the new variables zj obey the constraint (|18|) . Now taking into account that Q can be written as 

H^) = ( n \ = e''"^ , (22) 

it is not difficult to check that 

^^'K c d) = olt!!m,{s^}), (23) 

where 

, (24) 

KniV -U + v)Kn{r] - u' + v')Kn{u' - v)Kn{u - v') 

Further, using (jlSp and the reflection property d?]) for the elliptic gamma-function, one can re-write 
the ratio of the W-factors entering (I14p in the form 

_^ 2n 

W,,_,(b,a)W„._„(c,a)[w,._„(d,c)W.,._„(d,6)] = \{ T{tySk) (25) 

hk=l 



Next, substituting ([21]) into (j20|) one gets 

r. _ „-2(M'-i)')+2iCj r , _ 2{u-v)+2\bj ~. _ 2{n-(;'-»?)-2iaj ~ . _ 2{u' ^v-r{)-2\dj (c.r.-. 

Consider now the star weight (jl3p and make a change of variables Zj = e"^'^^ (note the minus sign 
in the exponent). Using the variables (I26p one obtains 



V-( ^^)=gIt'^{i^}ds^}) (27) 



where ^ is the same as in (j24p . The relations (|23p , (j25p and (|27p immediately imply that the star-star 
relation (|14p is equivalent to the Rains transformation formula (|19p for the elliptic hyper geometric 
integrals, obtained in [T9] . 
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